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Abstract 

We study the D — D mixing and rare D decays in the Littlest Higgs model. As the new weak 
singlet quark with the electric charge of 2/3 is introduced to cancel the quadratic divergence induced 
by the top-quark, the standard unitary 3x3 Cabibbo-Kobayashi-Maskawa matrix is extended to 
a non-unitary 4x3 matrix in the quark charged currents and .^-mediated flavor changing neutral 
currents are generated at tree level. In this model, we show that the D — D mixing parameter can be 
as large as the current experimental value and the decay branching ratio (BR) of D — > XyPj is small 
but its direct CP asymmetry could be 0(10%). In addition, we find that the BRs of D — > X u £ + £~ , 
D — > X u uv and D — > u + u~ could be enhanced to be O(10 -9 ), O(10" 8 ) and O(10~ 9 ), respectively. 
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I. INTRODUCTION 

As the observation of the B s — B s mixing in 2006 by CDF [jj , all neutral pseudoscalar- 
antipseudoscalar oscillations (P — P) in the down type quark systems have been seen. In 
the standard model (SM), the most impressive features of flavor physics are the Glashow- 
Iliopoulos-Maiani (GIM) mechanism Q and the large top quark mass. The former results 
in the cancellation between the lowest order short-distance (SD) contributions of the first 
two generations to the mass difference ArriK in the K° system, while the latter makes Atob ? 
(q = d, s) in the B q systems dominated by the SD effects [sj]. In addition, these features also 
lead to sizable flavor changing neutral currents (FCNCs) from box and penguin diagrams, 



which contribute to the rare decays, such as K — > ixvv and B — > 



t is known that 



3, 9. 



However, 



these processes could be good candidates to probe new physics effects 
it is clear that the new physics signals deviated from the SM predications for the P — P 
mixings and rare FCNC decays have to wait for precision measurements on these processes. 

Unlike K and B q systems, the SD contributions to charmed-meson FCNC processes, 
such as the D — D mixing j3] and the decays of c — > ui + £~ and D — > £ + £~ are highly 
suppressed due to the stronger GIM mechanism and weaker heavy quark mass enhancements 
in the loops. On the other hand, it is often claimed that the long-distance (LD) effect for 
the D — D mixing should be the dominant contribution in the SM. Nevertheless, because 



the nonperturbative hadronic effects are hard to contro 

US 



10 



11 



12|. Recently, BABAR 



the result is still inconclusive 



1. 



13|] and BELLE [lj, ll5| collaborations have reported the 



evidence for the D — D mixing with 

x' 2 = (-0.22 ± 0.30 ± 0.21) x 10" 3 , 

y' = (9.7 ±4.4 ±3.1) x 10~ 3 , (1) 

and 

x = = (0.80 ±0.29 ±0.17)%, 

r» 

A P 

y = — - = (0.33 ±0.24 ±0.15)%, 
y CP = (1.31 ±0.32 ±0.25)%, (2) 

respectively, where x' — x cos 5 + y sin 5 and y' = — x sin 5 + y cos 5 with the assumption 
of CP conservation and 5 being the relative strong phase between the amplitudes for the 



doubly- Cabbibo-suppressed D — > K + ti~ and Cabbibo-favored D — > K~n + decays [16|, Il7 ] 
and yep = t(D — > K~tt + )/t(D — ► K + K~) — 1. Moreover, no evidence for CP violation is 
found. The combined results of Eqs. (JTJ) and (J2]) at the 68% C.L. are [l8] 



x 

y 
5 



(5.5 ±2.2) x KT 3 
(5.4 ± 2.0) x 10~ 3 
(-38 ±46)°. 



(3) 



Note that the u ppe r bound of x < 0.015 at 95%C.L. can be extracted from the BELLE 



data in Eq. 
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151 ]. The evidences of the mixing parameters by BABAR and BELLE 



collaborations reveal that the era of the rare charmed physics has arrived. The results in 



Eq. (J3j) can not only test the SU(3) breaking effects 



examine new physics beyond the SM 17 
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2l|. 



or the D — D mixing but also 



It is known that a straightforward way to enhance the rare D processes is to include some 
new heavy quarks within the framework of the SM. For instance, if a new heavy quark with 
the electric charge of —1/3 is introduced, it could affect the D system since the extra down 
type quark violates the GIM mechanism. However, the constraint on this heavy quark is 
quite strong as it could also lead to FCNCs for the down type quark sector at tree level, 
which are strictly limited by the well measured rare K and B decays, such as Kl — > /U + (i~ 



and B — »• X s 7 22J. On the other hand, if the charge of the new heavy quark is 2/3, it could 
generate interesting tree level FCNCs for the up type quark sector, for which the constraints 
are much weaker. In this paper, we will study D physics based on a new weak singlet upper 
quark. 



It has been known that in the framework of the Littlest Higgs model 23], there exists 
a new SU(2)l singlet vector-like up quark 24j, hereafter denoted by T. Since the number 
of down type quarks is the same as that in the SM, the standard unitary 3x3 Cabibbo- 
Kobayashi-Maskawa (CKM) matrix is extended to a non-unitary 4x3 matrix in the charged 
currents. Moreover, Z-mediated FCNCs for the up quark sector are generated at tree level. 
In Ref. {25)] . it has been shown that the contributions of this new quark to the rare D 



processes are small and cannot reach the sensitivities of future experiments 



25 



26j. In this 



paper, we will demonstrate that by adopting some plausible scenario, the effects could not 



only generate a large D — D osci 
by BESIII for the rare D decays 



lation but also marginally reach the sensitivity proposed 



271 ] . We note that the implication of the new data on the 



D — D mixing in the Littlest Higgs model with T-parity has been recently studied in Ref. 
21]- 

The paper is organized as follows. In Sec. [Ill we investigate that when a gauge singlet T- 
quark is introduced in the Littlest Higgs model, how the non-unitary matrix for the charged 
current and the tree level Z-mediated FCNC are formed. By using the leading perturbation, 
the mixing matrix elements related to the new parameters in the Littlest Higgs model are 
derived. In addition, we study how to get the small mixing matrix element for V u ^b, which 
describes the b — > u(c) decays. In Sec. IHIt we discuss the implications of the non-unitarity 
on the D — D mixing and rare D decays by presenting some numerical analysis. Finally, we 
summarize our results in Sec. (TV] 



II. NON-UNITARY MIXING MATRIX IN THE LITTLEST MODEL 



To study the new flavor changing effects in the Littlest Higgs model, we start by writing 
the Yukawa interactions for the up quarks to be 



24 



25| 



£>Y — -labftijktxyXai'Ejx'F'kyU'b + l§fTT C + fl.C. , 

where xl — (di,ui,Q), x\ — ( s 2,c 2 ,0), x\ — (b 3 ,t 3 ,T), u' b is the weak singlet and E 



(4) 



e^Ene™^ with 



E 



/ 



L 2x2 



\ 



1 



2x2 



/ 



n 



h/V2 h*/V2 
h T /V2 



\ 



(5) 



The scale / denotes the global symmetry spontaneously breaking scale, which, as usual, 
could be around 1 TeV. Consequently, the 4x4 up-quark mass matrix is given by [25| 



Mi 



u 



ilijV 







3x3 



(6) 



l 3 sf | l f 

We remark that the quadratic divergences for the Higgs mass from one-loop diagrams in- 
volving t and T get exactly cancelled as shown in Ref. 28|. Moreover, for other quarks 
other than the top quark, the one-loop quadratic divergent contributions do not necessitate 
fine-tuning the Higgs potential as the cutoff is around 10 TeV for / ~ 1 TeV due to the small 



corresponding; Yakawa couplings. That is why there is no need to introduce extra singlet 



states T 



2n 



29) 



To obtain the quark mass hierarchy of m t ^> m c ^> m u , we can choose a basis such that 



the up-quark mass matrix is 



30] 




Mu = ( 7) 



where rhuij = 5ijliv/y2 = rrii is diagonal matrix and h = (hi, h 2 , h 3 ). The hi is related 
to Z 33 by hi = V^ 3 fl / 33 and hh^ = | Z33 1 2 , in which V Ur is the unitary transformation for the 
right-handed up quarks. We note that are not the physical masses and in principle their 
magnitudes could be as large as the weak scale. In order to preserve the hierarchy in the 
quark masses, one expects that m 3 > m 2 > mi. Furthermore, in terms of this basis, the 
charged and neutral currents, defined by 

C c = -i-J~W^ - -J- J-Wt» + h.c. , 

y/2 V2tan6 " H 



C N = (J3 - sin 2 6 w J» m ) Z, + t^-AZu, + h.c. , (8) 

cos oy/ tan u 

are expressed by 



^ = U Llll V°a v D L , 

J3 = \uL^V a v V^U L - l -D L YD L , (9) 
respectively, where U T = (ui, c 2 , t 3 , T), D T = (d, s, b), a v = diag(l, 1, 1,0) and 

V v 

with Vq KM Vq km = 13x3- The null entry in ay denotes the new T-quark being a weak 
singlet; and without the new T-quark, V^ KM is just the CKM matrix. Since the down quark 
sector is the same as that in the SM, we have set the unitary transformation U Dl to be an 
identity matrix. 

For getting the physical eigenstates, the mass matrix in Eq. (J7|) can be diagonalized by 
unitary matrices V Ul - r so that we have 

M diag M tdiag = V U LMuM ^yUl (U ) 



and 

MuM r = u u . 12 

\hrhy (i/ 33 i 2 + i/oi 2 )/ 2 ; 

Since ( | ^33 1 2 + \lo\ 2 )f 2 is much larger than other elements, we can take the leading order of 
the perturbation in hiirii/ f (i = 1,2,3). According to Eq. (ITT]) , the leading expansion is 
given by 

M diag M Jiiag = V Ul MuM \ uV U{ „ (1 + A £ )M^Mj(l - A L ) . (13) 

By looking at the off-diagonal terms (M^ iag M^ diag )j4(4j), we can easily get 

A ~ A - himf fl/fl 

with % 7^ 4. From the diagonal entries, if we set the light quark masses to be m u m c 0, 
we obtain 

« m 2 u . « m) + 2A Lj4 (M u Ml 7 ) Aj , 

" "2ft"? (15) 
with j = 1, 2. To be consistent with Eq. (114)) . at the leading expansion the relation 

2h) = (|fe| 2 + |/o| 2 ) (16) 

should be satisfied. We emphasize that the choice of Eq. (ITB]) is somewhat fine-tuned in 
order to have Eqs. (114j) and (115j) simultaneously. Since the top-quark is much heavier than 
other ordinary quarks, we have 2h\ ~ (1 — / ?n§) ( | Z33 ] 2 + |/o| 2 ) if / > m 3 > m <- Similarly, 
one obtains the flavor mixing effects for i ^ j ' ^ 4 to be 



AUJ ™\ - mj ((|/ 33 | 2 + |/ | 2 )/ 2 - m|) ((|/ 33 | 2 + |/o| 2 )/ 2 - ml) ' 
After diagonalization, the currents become 

J- = U Llfl V UL V°a v D L = U Lltl V UL V°D L , 

= U L ^V UL V°a v V 0] V ]UL U L = U L YV UL V°V oj V WL U L , (18) 

where U T = (u, c, t, T), 

yO = yO av = I (V8KM)S X3 





and diag(V A °V A °^) = ay- Since the 4-th component of ay is different from the first 3 ones, it 
is obvious that the matrix V = V Ul V° associated with the charged current does not satisfy 
unitarity. In addition, V UL ayV u L , which is associated with the neutral current, is not the 
identity matrix. As a result, Z-mediated FCNCs at tree level are induced. According to 
Eq. (ITS]) , we see that 

W* = V UL a v V u L (20) 

which is just the same as the effects of the Z-mediated FCNCs. Due to V being a non-unitary 
matrix, one finds 

{VV% = V i4 V* A . (21) 

Consequently, the interesting phenomena arising from non-unitary matrix elements are al- 
ways related to VnV* A = A I/i4 Aj 4 . We note that as we do not particularly address CP 
problem, in most cases, we set the parameters to be real numbers. 

It has been known that enormous data give strict bounds on the flavor changing effects. 
In particular, the pattern describing the charged current has been fixed quite well. Any 
new parametrization should respect these constraints. It should be interesting to see the 
relationship with and without the new vector-like T-quark. From Eq. (Tl8l) . we know that 
the new flavor mixing matrix for the charged current is given by V = V Ul V°. At the leading 
order perturbation, one gets 

V = V Ul V°& (1 + A L )V° = V° + A L V° . (22) 

If V t ° b ~ 1 is taken, one finds that V ub « V° b + A L iz and V cb « V° b + A L23 . In terms of Eq. ([TTJ) 
and hi — /i 2 ~ h%, the relations A^ 13 w —mi/m 3 and Ax 23 m — m 2 /m 3 are obtained. Hence, 
in our approach, we have 

, V^Vg- — - (23) 
m 2 m 3 m 3 

From these results, it is clear that when the T-quark decouples from ordinary quarks, 
V us -> V° s , V ub V° b and V cb -> V^, while m 1 /m 2 -> m u /m c , m x jm z m u /m t and 
m 2 /m 3 — > m c /m t , respectively. According to the observations in the decays of b — > uli>i and 
b — > c£ue, the corresponding values have been determined to be |T4&| = 3.96 ± 0.09 x 10~ 3 



and \V cb \ = 42.21±g;|g x 10" 3 , respectively 



161 ] . Since and m ; are free parameters, to 



satisfy the experimental limits with interesting phenomena in low energy physics, it is rea- 
sonable to set the orders of magnitude for mi/m^ and m^jvn^ (rai/7712) to be O(10 -2 ) and 
O(10 _1 ), respectively. Consequently, the non-unitary effects on the rare charmed meson 
decays governed by VuV^ could be as large as A 14 A 2 4 ~ O(10 -4 ), which could be one order 



of magnitude larger than those in Ref. 251 ] . 



III. D-D MIXING AND RARE D DECAYS 



A. D — D mixing 



It is well known that the GIM mechanism has played an important role in the K — K 
oscillation in the SM. In addition, due to the top-quark in the box and penguin diagrams, 



B q — B q mixings are dominated by the SD effects, which are consistent with the data 16] 



On the contrary, for the D — D mixing the GIM cancellation further suppresses the mixing 
effect to be Am^ ~ 0{m A s /m 2 v m 2 c ) |7| and the bottom quark contribution actually is a 
subleading effect due to the suppression of {V u bV* b ) 2 . In the SM, the SD contribution to the 
mixing parameter is O(10 -7 ) 34(. However, the LD contribution to the mixing is believed to 
be dominant. Due to the nonperturbative hadronic effects, the result is still uncertain with 

n _ n n n 

the prediction on the mixing parameter ranging from O(10" 3 ) |9| to O(10" 2 ) [Ifl, llll, ll2|. 
Nonetheless, the mixing parameters shown in Eq. ([3]) could arise from the LD contribution. 
Thus, it is important to have a better understanding of the LD effect. On the other hand, it 
is also possible that the mixings in Eq. ([3]) could result from new physics. In the following, 
we will concentrate on the Littlest Higgs model. 

In the quark sector of the Littlest Higgs model due to the introduction of a new weak 
singlet, a direct impact on the low energy physics is the FCNCs at tree level. According to 
Eq. (1TB]) . the most attractive process with |AC| = 2 via the Z-mediated c—u — Z interaction, 
illustrated in Fig. [TJ, is given by 

W(|AC|=2) = -^(V u V 2 \) 2 u^P L cuYP L c, 

F (VuV* 4 ) 2 u ltl P L cuYP L c (24) 



V2 

In terms of the hadronic matrix element, defined by 



(D\(uc) v _ A (uc) v _ A \D) = -B D f 2 D m 2 D , (25) 
8 




FIG. 1: Z-mediated flavor diagram with |AC| = 2. 



the mass difference for the D meson is 25] 

V2 



Am 



i) 



G F f 2 D m D B D \(V u V 2 \ 



* \2| 



(26) 



If we assume no cancellation between new physics and SM contributions, by taking = 
1/T D = 6.232 x 10 11 GeV" 1 , = 200 MeV [3ll, Isj and m D = 1.86 GeV and using 

Eq. ([3]), we obtain 



Co = \V U V* 4 \ = |A L14 A L24 | = (1.47 ±0.29) x 10" 



(27) 



which is in the desirable range. In other words, the result in Eq. (I27p demonstrates that the 
non-unitarity in the Littlest Higgs model could enhance the D — D mixing at the observed 
level. We note that the limit of x < 0.015 (95%C.L.) leads to 



Co < 2.5 x 10" 4 . 



(28) 



In addition, we note that cancellation between the LD effect in the SM and the SD one from 
new physics could happen. In this case, the values in Eqs. (!27|) and (|28|) could be relaxed. 



B. 



D X u -y decay 



In the SM, the D-meson FCNC related processes are all suppressed since the internal 
fermions in the loops are all much lighter than my/. For the decay of D — >• X M 7, without 
QCD corrections, the branching, ratio is O(10~ 17 ); and it becomes O(10~ 12 ) when one-loop 
QCD corrections are included 8J. However, it is found that the two-loop QCD corrections 



can boost the BR to be as large as 3.5 x 10' 



. It should be interesting to see how large 



the non-unitarity effect on c — > w~f is in the Littlest Higgs model. 



9 



To study the radiative decay of c — > u^y, we write the effective Lagrangian to be 

G Frr e 



£ c ^u~f = -—^VusV^Cy—rricUa^PRcF 



(29) 



33] 



where C 7 = C 7 5M + C? p and Cf M « -(0.007 + z'0.02) = 0.021e i,5s with 5 S = 70.7° 
being the strong phase induced by the two-loop QCD corrections. In the extension of the 
SM by including a weak singlet particle, the flavor mixing matrix in the charged current 
is not unitary and the Z-mediated FCNC at tree level is generated as well. For c — > uj, 
besides the QED-penguin diagrams induced by the IV-boson displayed in Figs. [2ti and [2b, 
the Z-mediated QED-penguin one in Fig. will also give contributions. We note that the 

d,s,b W Z 




FIG. 2: Flavor diagrams for c — ► wy. 



contributions from Wh and Zh can be ignored as m 2 v /m 2 VH and m 2 z /m 2 ZH are much less 
than one. 

At the first sight, due to the light quarks in the loops, the contributions from Figs. ^ 
and [2b could be negligible. However, due to the non- unitarity of (VV^) UC = V 14 V 2 \ ^ 0, 
even in the limits of m^ Sj 6 — > 0, the contributions from the mass independent terms do not 
vanish anymore and can be sizable. In terms of unitary gauge 22J, we obtain 

ll(Wt) 12 11V U V& 



Furthermore, if we set m u 



36 V US V* 



36 V US V* 



(30) 



0, the contributions from Fig. &p are given by 



fc + fu 



i \ (i 

7 1 1 1 



e u sm 2 e w ) [4^(0) -66(0) + 26(0)] 



2 { \2 

+e u sm 2 6 w [4^(0) -46(0)]} , 

\e u [26(2/t)-36(2/t)+6(2/t)] 



where the functions £ n (x) are defined by 



z n+1 dz 
o l + (x-l)z 



(31) 



(32) 



10 



and yT = rriT/mz- Numerically, the total contribution in Fig. [2] is 

0.53 



w 

7 



v*v 

v cs y us 



v u v* 4 



(33) 



If we regard V14K24 as an unknown complex parameter, i.e. V^V^ = (oe 10 with 9 being 
the CP violating phase, one can study the decay BR and direct CP asymmetry (CPA) of 
D — > X u j defined by 



BR(D -> X ul ) 
A C p 





c 7 


2 


71 


v cd 


2 



BR(D -> X d eu e 



BR(c -> tZ7) - BR(c -> 1*7) 



(34) 



BR(c -> S7) + BR(c -> «7) 

as functions of ^ an d ^- 111 Fig- El the BR and CPA as functions of Co are presented, 
where the solid, dotted, dashed and dash-dotted lines represent the CP violating phase at 
9 = 0, 45°, 90° and 135°, respectively. From these results, it is interesting to see that 

3.7 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 15 r 




"0 1 2 3 4 5 6 
(a) S ° 



1 2 3 4 5 6 

(b) 



FIG. 3: BR (in units of 10 -8 ) and CPA (in units of 10~ 2 ) for D — ► Xuj as functions of £0, where 
the solid, dotted, dashed and dash-dotted lines represent the CP violating phase at 6 = 0, 45°, 90° 
and 135°, respectively. 



BR(D — > Xu'j) is insensitive to the new physics effects, whereas the direct CPA could be as 
large as O(10%). Explicitly, if we take 9 = 90° and Co = 1.5 X 10" 4 , the CPA is about 3%. 
Note that this CPA vanishes in the SM. 



11 



C. D -> X u ££ and D° -» £+£~ decays 

Because the current experimental measurements in K and B q decays are all consistent 
with the SM predictions, it is inevitable that if we want to observe any deviations from the 
SM, we have to wait for precision measurements for K and B q . SuperB factories or LHCb 
could provide a hope. However, the situation in D physics is straightforward. As stated 
before, unlike K and B q systems, due to no heavy quark enhancement in the D system, the 
rare D-meson decays, such as D — > X u ££ (£ = e, /x, u), are always suppressed. Even by 
considering the long-distance effects, the related decays, such as D — > fi + n~ and D — > X u i/P, 



get small corrections to the SD predictions on the BRs 35j. Therefore, these rare decays 



definitely could be good candidates to probe the new physics effects. Since the values in the 
SM are hardly reachable at D factories [27], if any exotic event is found, it must be a strong 
evidence for new physics. In the following analysis, we are going to discuss the implication 
of the Littlest Higgs model on the rare D decays involving di-leptons. 

To study these decays, we first write the effective Hamiltonian for c — > u£ + £~ {£ = e, /i) 

as 

H( C - u£ + r) = -^^V C * S V US [Cf0 9 + C 7 7 + C 10 O l0 ] , (35) 

V 2tt 

2 TJl _ 

o 7 = — ^mo^p R chn, 

9 = U^P L drfl, 

O10 = u^P L ch^ 5 l, (36) 
where the effective Wilson coefficients are given by 

Cf = —^^4 + (h(z„ s) - h(z d , s)) (C 2 (m c ) + S^K)) , 

10 ~ ~a \W Al { } 

u em ' cs' us 

with s = q 2 /ml, z\ = mi/m c , Cy = — 1/2 + 2 sin 2 9w, c e A = —1/2 and 



4 4 2 1 1— 
h(z,s) = — Inz + — + -x (2 + x)v|l 



9 27 9 9 



In 



x 



0=x-l 



(38) 



i it, for x = 4z 2 / s < 1 , 
2 arctan^=, for x = 4z 2 /s > 1 

Vi- 1' ' 

Here, we have neglected the small contributions from the penguin and box diagrams. We 



note that in the SM, the SD contributions are mainly from the term with h(z, s), induced 

12 



by the insertion of O2 = ul^^l^lI^cl and mixing with O9 at one-loop level [35|, l36 |. 
We note that the resonant decays of D — > X U V — > X u £ + £~ (V = <p, p, uj) would have 
large corrections to c — > u£ + £~ at the resonant regions. However, in this paper we do 
not discuss these contributions as we only concentrate on the SD contributions. Moreover, 
these resonance contributions can be removed by imposing proper cuts in the phase space 
in dedicated searches. 

From Eq. fl35l) . the decay rate for D — > X u £ + £~ as a function of the invariant mass 
s = q 2 /m 2 can be found to be 

dT G 2 F mla 2 cm 2 _ 2 

ds 768tt5 ]VusVcsl [i S) K[S) > 

R(s) = (\Cf I 2 + ICiol 2 ) (1 + 2s) + 12Re(C* 7 Cf ) + 4 M + |C 7 | 2 . (39) 

In addition, by utilizing the lepton angular distribution, we can also study the forward- 
backward asymmetry (FBA), given by 

dApB f-i d cos GdY /dsd cos 6 sgn(cos 6) 



ds d cos 8dT/ dsd cos 8 



-3—rx Re 
R(s) 



Cf + -C 7 ) C* 



(40) 



where 6 1 is the angle of £ + related to the momentum of the D meson in the £ + £~ invariant 
mass frame. Since C\q is small in the SM, Aps is negligible. With m c = 1.4 GeV and the 
mixing parameter in Eq. ff27j) . we get 



BR(D — > X u e + e~) = (4.18 ± 0.91) x 10" 10 , 

BR(D -> X u fx + fi~) = (2.51 ± 0.86) x 10~ 10 , (41) 

comparing with the SM predictions of BR(D — > X M e + e")5M = 2.1 x 10 -10 and BR(D — > 
X u n + n~) S M = 0.5 x 10^ 10 , respectively. Clearly, if some cancellation occurs between new 
physics and SM contributions in the D — D mixing, a larger value of Co can be allowed. In 
Fig. HI we show the tendency of the decay as a function of (0, where the negative horizontal 
values correspond to -Co- In addition, we present the differential decay BR [FBA] of 
D — > X u e + e~ as a function of s = q 2 /m 2 c in Fig. \5$i [b], where the thick solid, dotted and 
dashed lines correspond to Co = 1-5, 3.0 and 5.0, while the thin ones denote the cases for 
—Co except Co = for the thin solid line in Fig. [5^. From Fig. [5b, we see that the FBA 
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t_l_i 





4 



FIG. 4: BR(in units of 10 9 ) for D — ► X u e + e as a function of Co- 

0.2 




0.2 0.4 0.6 0.8 

s 

(a) 



0.2 0.4 0.6 0.8 

s 

(b) 



FIG. 5: (a) [(b)] Differential BR (in units of 1CT 9 ) [FBA] for D — > X u e + e~ as a function of s, where 
the thick solid, dotted and dashed lines correspond to Co = 1-5, 3.0 and 5.0, while the thin ones 
denote the cases for — Co except Co = for the thin solid line in (a). 

is only at percent level. In the Littlest Higgs model, this is because the Z coupling to the 
charged lepton Cy = — 1/2 + 2 sin 2 6^ appearing in Cg S is much smaller than one. This is 
quite different from that in b — > s£ + £~ where the dominant effect in the SM for the FBA is 
from the box and QED-penguin diagrams. 

Next, we discuss the decay of D — > X u vv. In the SM, the BR for D — > X u vv is estimated 
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to be O(10 16 ) — 0(10 15 ) 35j, which is vanishing small. In the Littlest Higgs model, by 
taking C 7 = 0, Cf = -C 10 = -n{VV^) 14 /{a cm V us V c * s ), the effective Hamiltonian in Eq. ([35]) 
can be directly applied to c — > uvv. The decay rate for D — > X u vv as a function of s = q 2 /m 2 
can be obtained as 



dT 

ds 



s) 2 (l + 2s) (27r 2 |(Wt )l2 |2) 



(42) 



768tt 5 v 

where the factor of 3 stands for the neutrino species. With Co = 1.5 x 10~ 4 , we get BR(D — > 
XuUu) = 1.31 x 10 -9 . However, if we relax the constraint on VuV^, the BR as a function 
of Co is shown in Fig. [6^- For a larger value of Co, BR(D — > X u vv) could be as large as 

orio- 8 ). 




FIG. 6: (a) BR (in units of 10" 8 ) for D -> X u vD and (b) BR (in units of 10~ 9 ) for D -> /i+fT. 

Finally, we study the decays of D — > . It has been known that, in the SM, the SD 
contributions to D — > are O(10~ 18 ), while the LD ones are O(10 -13 ) 35J. It is clear 

that any signal to be observed at the sensitivity of the proposed detector, such as BESIII, 
will indicate new physics effects. Since the effective interactions at quark level are the same 
as those in Eq. ( l35i) . one finds that 



BR(D 



G - 2 < 2 -A Am h,v u v 2 M 2 



16n'- 



;T D m D m e f D Jl - 



rn 



(43) 



D 



Here we have used equation of motion for the charged lepton so that £$ D £ = 0. We also 
note that operators 0-j$ make no contributions. With (V^V^I = Co = 1-5 x 10~ 4 , the 
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predicted BR for D — > fi + fi is 1.17 x 10 10 . In Fig. [6b, we present the BR as a function of 
Co- We see that BR(D — > in the Littlest Higgs model could be as large as O(10~ 9 ). 



IV. CONCLUSIONS 

We have studied the D — D mixing and rare D decays in the Littlest Higgs model. In 
the model, as the new weak singlet vector- like quark T with the electric charge of 2/3 
is introduced to cancel the quadratic divergence induced by the top-quark, the standard 
unitary 3x3 CKM matrix is extended to a non-unitary 4x3 matrix in the quark charged 
currents and Z-mediated flavor changing neutral currents are generated at tree level. We 
have shown that the effects on |AC| = 2 and |AC| = 1 processes are all related to VmV^ in 



Eq. (ED. 



To avoid the scenario adopted by Ref. 25J, in which Iq ~ £33 ^> was assumed, we 



choose the basis such that the effective mass matrix for ui, C2 and £3 is diagonal, while 
the corresponding masses m 1; m 2 and m 3 are free parameters and can be as large as the 
weak scale v. Since the global symmetry breaking scale / is larger than v, the mixing 
matrix relating physical and unphysical states could be extracted by taking the leading 
perturbative expansion. Accordingly, by using the approximation of m u m m c m 0, the 
explicit expressions for Vu and V24 have been obtained. In terms of the data for V u b and 
Vcb, we have found that the natural value for £ = |Vx 4 V^ 4 | is O(10~ 4 ), which agrees with the 
observed parameter in the D — D mixing but it is one order of magnitude larger than that 



in Ref. 



25|. 



For the rare D decays, due to the non-unitarity effects in the model, BR(D — > X u £ + £ ), 
BR(D -> X u vv) and BR(D -> ^V - ) could be enhanced to be O(10~ 9 ), O(10" 8 ) and 
O(10~ 9 ), respectively, which could marginally reach the sensitivity proposed by BESIII j^. 
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